ICCUB-11-002 
UB-ECM-PF 10/44 



OVERVIEW OF CHARMONIUM DECAYS AND PRODUCTION 
FROM NON-RELATIVISTIC QCD 



JOAN SOTO 

Departament d'Estructura i Constituents de la Materia and Institut de Ciencies del Cosmos, 
Universitat de Barcelona, Marti i Franques 1 
Barcelona, Catalonia, Spain 

joan.soto@ub.edu 



I briefly review Non-Relativistic QCD and related efli'ective theories, and discuss appli- 
cations to heavy quarkonium decay, and production in electron-positron colliders. 

Keywords: Heavy Quarkonium; Non-Relativistic QCD. 

PACS numbers: 12.38.-t, 14.40.Pq, 12. 39. St, 13.20.Gd, 13.25. Gv 



1. Introduction 

Heavy quarkonia are mesons made out of a heavy quark and a heavy antiquark 
whose masses (toq) are larger than Kqcd, the typical hadronic scale^. These in- 
clude bottomonia (66), charmonia (cc), Be systems (6c and c6) and would-be to- 
ponia (if). In the quarkonium rest frame the heavy quarks move slowly {v ^1, v 
being the typical heavy quark velocity in the center of mass frame), with a typical 
momentum toqu <C mq and binding energy raqv^ . Hence any study of heavy 
quarkonium faces a multiscale problem with the hierarchies mg 3> mqv 3> rnqv^ 
and mg ^ ^qcd- The use of effective field theories is extremely convenient in order 
to exploit these hierarchies. Based on the pioneering work of Caswell and Lepagj2l^ 
a systematic approach to study these systems from QCD has been developed which 
is generically known as Non-Relativistic QCD (NRQCDj^l. The NRQCD formalism 
for spectroscopy and inclusive decays to light particles is very well understood (this 
is also so for electromagnetic threshold production) so that NRQCD results may be 
considered QCD results up to a given order in the expansion parameteri^ (usually 
as(mg) and 1/mg). The NRQCD formalism for production is more controversial. 
I will summarize here the main features of NRQCD (and related effective theories) 
and review recent theoretical progress in decays, and production at electron-positron 
colliders. I refer to the contribution of Geoff Bodwin for an overview on production 
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at electron-proton and hadron colliders, and on the status of factorization proofi^. 

2. Non-Relativistic QCD 

The hierarchy of scales exploited in NRQCD is mg ^ mQV,mQv'^ , Aqcd- The 
part of the NRQCD Lagrangian bilinear on the heavy quark fields coincides with 
the one of Heavy Quark Effective Theory (HQET)^, but the size assigned to each 
term differs from it. It also contains analogous bilinear terms for the antiquark, 
and four fermion operatori^, which are classified as color singlet or color octet. The 
short distance matching coefficients of these operators have imaginary parts. The 
fact that NRQCD is equivalent to QCD at any desired order in l/mg and asimg) 
makes the lack of unitarity innocuous. In fact, it is turned into an advantage: it 
facilitates the calculation of inclusive decay rates to light particles. 

The inclusive decay widths of heavy quarkonium states to light particles are 
given in terms of sums of NRQCD matrix elements of both color singlet and color 
octet operators multiplied by the imaginary part of short distance matching coef- 
ficients, which can be calculated in perturbation theory in Q!s(toq). Once a size in 
terms of mg and v is assigned to each of the matrix elements, a systematic expan- 
sion in powers of v and as{mQ) is obtained. Earlier QCD factorization formulas 
were missing the matrix elements of color octet operators. They are inconsistent 
because the color octet matrix elements are necessary to cancel the factorization 
scale dependence which arises in loop calculations of the short distance matching 
coefficientJ^. The matrix elements cannot be calculated in perturbation theory of 
Q!s(wq). The color singlet ones are related to wave functions at the origin but the 
color octet ones are not, in the general case. Sometimes only the effect of color o ctet 
operators in the renormalization group evolution equations is taken into accounlPISl. 
If niQV A.QCD they can be estimated using weak coupling techniques EHH! and 
if rnqv^ <IC Aqcd they can be related to wave functions at the origin plus a num- 
ber of universal non-perturbative parameter&li^. They can alwa ys b e calculated on 
the lattice-*^^. Alternatively they can also be extracted from datsll2l. The imaginary 
parts of the ma tchi ng coefficients of the dimension 6 and dimension 8 operators are 
known at 0{al^, those of dimension 10 (g and P waves) at (^(q;^ j l5il6 | 17 | 18 | 19 | 
those of dimension 10 {D wa ve) at 0(03 j^SEEl]^ ^^^^ -^j^g ^^^g ^y^e color singlet 
operator at 0{at)^^. 

3. Potential NRQCD 

Unlike HQET, the NRQCD Lagrangian does not enjoy a homogeneous counting, 
due to the fact that the scales mgu, toqu^ and Aqcd are still entangled. Potential 
NRQCD (pNRQCD) aims at disentangling these scales, and hence to facilitate the 
counting, by constructi ng a further effective theory in which energy scales larger 
mgu^ are integrated oull^. If Aqcd ^ "mqv'^ , then rnqv ^ Aqcd and the match- 
ing between NRQCD and pNRQCD can be carried out in perturbation theory in 
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as{mQv) . This is the so called weak coupling regime. If Aqcd ^ i^^Qv"^, the match- 
ing cannot be carried out in perturbation theory in as{mQv) anymore, but one can 
still exploit the hierarchy mg » mqv.KQCD ^ rnQv"^. This is the so called strong 
coupling regime. Since mg, v and Aqcd are not directly observable, given heavy 
quarkonium state it is not clear to which of the above regimes, if to anjQ , it must 
be assigned to. A test was proposecP^ using the photon spectra of radiative decays. 
Under this test, CLEO datJM suggests that T{2S) and T{3S) are in the strong 
coupling regime whereas T(IS') is not. An alternative test using leptonic decays 
was also carried out leading to the same conclusiorPZl. Even though experimental 
data is now available"^^ and it was suggested at some poinlP^, no such a test for 
charmonium states has been carried out yet. 



3.1. Weak Coupling Regime 

The pNRQCD Lagrangian in this regime can be written in terms of a color sin- 
glet and a color octet heavy quark-antiquark wave function fields that evolve ac- 
cording to the Hamiltonians hs and ho respectively, and interact with gluons of 
energy ^ mQV^ (ultrasoft). The potentials in those Hamiltonians, as well as the 
remaining matching coefficients, may be obtained by matching to NRQCD in per- 
turbation theory in as{mQv) and l/rng at any order of the mu ltipole expans^ 
(l/mgw). The static potential in hs is known up to three loop j^^ l ^^ l ^'^ l ^^ l ^^ l ^ ^ 
(the logarithmic contributions at four loops are also knowrPSl) ^ and in ho up to two 
loops . The renormalon singularities in the static potentials are also understood 
in some detaii^SI. The 1 /mq and 1 /wg terms in hs are known at two and one loop 
respectiveljSH. This Lagrangian has been used to carry our calculations at fixed 
order in a compl ete NN NLO (assuming Aqcd ^ ''^Q'^D expression for the 
spectrum is available | 40|35| ]y[ost remarkably resum mations of logarithms can also 
be carried out using renormalization group techniqueJSHHSl. Thus the hyperfine 
splitting and the ratio vector /pseudoscala r elect romagnetic decay width have been 
calculated at NNLL and NLL respectiveljEESl, xhe hyperfine splitting compares 
well to data for charmonium, but undershoots considerably the experimental value 
for bottomonium. A reorganization of the quantu m m echanical calculation has been 
recently proposed that may resolve this problerrPSI^ At least it brings the ratio of 
the vector/pseudoscalar electromagnetic decays mention above in good agreement 
with experiment for charmonium. 



3.2. Strong Coupling Regime 

The pNRQCD Lagrangian in this regime reduces to a heavy quark-antiquark wave 
function field interacting with a potentialB, which can be organized in powers of 

^States close or above the open flavor threshold are expected to belong neither the weak nor the 
strong coupling regimes. 

''We ignore for simplicity pseudo-Goldstone bosons. 
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l/mg. These potentials cannot be calculated in perturbation theory of as{mQv) 
anymore, but can be calculated in lattice simulation J^Zl. They must coincide at 
short distances with perturbative calculationJ^S! and at long distances they must be 
compatible with the effective string theory of The fact that Aqcd S> mqv 

can now be exploited to further factorize NRQCD decay matrix element s into wave 
functions at the origin and universal bound state independent parameter j-^-^l^^ l. This 
further factorization allows to put forward new model independent predictions, for 
instance the ratios of hadronic decay widths of P-wave states in bottomonium were 
predicted from charmonium dat A or the ratio of photon spectra in radiative 
decays of vector resonanceJ^H. 



4. Beyond Inclusive Decays 

NRQCD can in principle describe transitions between heavy quarkonium states, 
but no useful information has been extracted so far from it, beyond the fact that 
non-relativ istic had ronic effective theories implementing heavy quark and chiral 
sy mmet rie j ^ ^ I ^ ^ may be considered hadronic realizations of it. pNRQCD can 
also describe transitions between heavy quarkonium states below open flavor thresh- 
old, and has been applied to magnetic dipole transitioni^. Decays to heavy light 
meson pairs can in principle be described by NRQCD as well, but again no much 
information has been obtained from it so far. 

For exclusive decays and for certain kinematical end-points of semi-inclusive 
decays to light particles, NRQCD must be supplemented with collinear degrees of 
freedom. This can be done in the effective theory framework of SCEl^^. Exclusive 
radiative decays of heavy quarkonium in SCET have been addressed^, where re- 
sults an alog ous to those of traditional light cone factorization formulas have been 
obtained^ Recently, SCET, HQET, NRQCD and pNRQCD have been combined to 
obtain model independent formulas for bottomonium decays to DD meson pairJSH 
Concerning semi-inclusive decays, a remarkable success has been achieved by com- 
bining SCET, NRQCD and pNRQC D for th e photon spectrum of T(15') 7X, a 
challenge since the early days of Unfortunately, when this forma hsm is 

applied to radiative decays of J/ tJj the outcome suffers from large uncertaintieil^^^. 



5. Production 

In production processes, like in semi-inclusive and exclusive decays, gluons of en- 
ergies ^ niQ may exists in the initial or final states, and hence NRQCD must be 
supplemented with additional factorization formulas in order to disentangle them 
from the dynamics of heavy quarkonium. 



5.1. Electromagnetic threshold production 

This is the simplest and best understood production process because no gluons of 
energies ~ toq exist in the initial and final states. In the weak coupling regime. 
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which is relevant for a precise measurement of topquark mass in the future ILC, 
the cross-section at NNLO is known for some time^ and the log resummation at 
NLL is also available™. The calculations at NNNLO are almost complet d ^ 
and partial results also exist for the log resummation at 

NNLlJMIMl. 

At this level 

of precision electroweak effects must also be taken into account in the tt casJ^. 



5.2. Inclusive production 

A factorization formula for the inclusive production of heavy quarkonium was put 
forward in the framework of which was assumed to hold provided that 

the transverse momentum was larger or of the order of the heavy quark mass. 
The formula contains a partonic level cross section, which can be calculated in per- 
turbation theory in as{mQ), in which the heavy quark pair may be produced in 
color singlet or a color octet state, and long distance production matrix elements 
which encode the evolution of the pair to the actual physical state, which cannot be 
calculated in perturbation theory in q;s(toq), and are usually extracted from data. 
The production matrix elements were assigned sizes according to the NRQCD veloc- 
ity scaling rules- , which correspond to the weak coupling regime. This factorizat ion 



formalism is receiving a closer look in the framework of fragmentation functions^ 
(px ^ Tig), and has been proved to be correct at NNLO in as(TOg), provided a 



slight redefinition of the matrix elements is carried oul)^. 

The measurement of a{J/ip + cc + X)/a{J/->p + X) at Belle in 200^0! -^as in 
clear conflict with LO NRQCD resulti^, and has triggered a number higher order 
calculations. For a(J/iJ;+cc+X), a NLO calculation of the color singlet contribution 
in asimq) g ives a large correctiorP2ll3| xhis is also the case for the NLO color octet 
contributiorPSl^. Relativistic correctionJ^ and two-photon contributioni^ have 
also been calculated. For a {J/il} + X(non— cc)), the LO color octelP^and the NLO 
color single both in (juq and contributions have been calculated. When 
all these results are put together there is no obvious discrepancy with experiment 
anymore, but the theoretical uncertainties are still large. 

Near certain kinematical end-points NRQCD production processes must be sup- 
plemented with collinear degrees of freedom, in analogy t o se mi-inclusive decays 



discussed above. This issue has been addressed using SCEtSSI^ 



5.3. Exclusive production 

The basic ideas of NRQCD factorization have also been applied to exclusive pro- 
cesses, mostly af ter the surprisingly large double charmonium cross-section first 
measured at BellJ^. Recently, factorization proofs have been put forward for this 
process as well as for exclusive production in i?-decayJ^. Although Belle and Babar 
results for J/ip+rjc production cross section do not quite agree with each other, b oth 
numbers are roughly an order of magnitude larger than LO NRQC D predictionJ^. 
A NLO calculation in a(rnQ) gives a very large correctiorPSEHI which, together 
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with a number of relativistic correctionJ^, brings th e theo retical calculation in 
agreement with both experimental results within errori^^SISSj 
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